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ABSTRACT 

We present the first part of the observations made for the Continuum Halos in Nearby Galax¬ 
ies, an EVLA Survey (CHANG-ES) project. The aim of the CHANG-ES project is to study 
and characterize the nature of radio halos, their prevalence, as well as their magnetic fields and 
the cosmic rays illuminating these fields. This paper reports observations with the compact 
D-configuration of the Karl G. Jansky Very Large Array (VLA) for the sample of 35 nearby 
edge-on galaxies of CHANG-ES. With the new wide bandwidth capabilities of the VLA, an un¬ 
precedented sensitivity was achieved for all polarization products. The beam resolution is an 
average of 9.6'^ and 36" with noise levels reaching approximately 6 and 30 /rJy/beam for C and 
L bands respectively (robust weighting). 

We present intensity maps in these two frequency bands (C and L), with different weightings, 
as well as spectral index maps, polarization maps and new measurements of star formation rates. 
The data products described herein are available to the public in CHANG-ES Data Release 
available at http://www.queensu.ca/changes, 

We also present evidence of a trend between galaxies with larger halos having higher SFR 
surface density, and show, for the first time, a radio continuum image of the median galaxy, 
taking advantage of the collective signal-to-noise of 30 of our galaxies. This image shows clearly 
that a ’typical’ spiral galaxy is surrounded by a halo of magnetic fields and cosmic rays. 

Subject headings: radio continuum: galaxies — galaxies: surveys, ... 
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1. Introduction 


This is the fourth paper in the series ’’Contin¬ 
uum Halos in Nearby Galaxies - an EVLA Survey” 
(CHANG-ES, Irwin et al. 20121. The overall aims 
of the GHANG-ES project are to investigate the 
occurrence, origin and nature of radio halos, to 
probe the disk-halo interface, and to investigate 
in-disk emission. We seek to understand connec¬ 
tions between radio halos and the host disk and 
its environment, and to investigate the magnetic 
fields in these galaxies and their halos. This cov¬ 
ers a wide variety of science, including cosmic ray 
transport and wind speeds, the nature and origin 
of galactic magnetic fields, the far-infrared radi¬ 
ation - radio continuum correlation, cosmic rays 
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and high-energy modelling, disk star formation 
rates and presence or absence of active galactic 
nuclei (AGN). Note that in this context, we use 
the word ’’halo” to refer to gas, dust, cosmic rays 
and the magnetic field above and below the galaxy 
disk, not to be confused with stellar or dark mat¬ 
ter halos. Specifically, we call emission on larger 
scales, i.e. scale height z > 1 kpc, halo emission, 
while the disk-halo interface is at 0.2 < z < 1 kpc 
(as defined in Irwin et al. [2012 ). 

Magnetic fields in the halos of edge-on galaxies, 
and their commonly X-shaped behaviour (mag¬ 
netic field lines outside the projected galaxy’s disk 
are bending away from the disk with increasing 
vertical components) are discussed in the review 
by Krause (2009) (also see references therein). 
Their intensities and degree of uniformity can 
be estimated from the total and polarized syn¬ 
chrotron emission, but many uncertainties per¬ 
taining to their structures and origins in external 
galaxies remain (Haverkorn fc Heesen||2012 1. 


With CHANG-ES, we have observed thirty-five 
nearby edge-on galaxies in the radio continuum in 
L and C bands (centred at approximately 1.5 and 
6 GHz, respectively), in three array configurations 
(B, C, D; in B-configuration only L-band was ob¬ 
served) of the Karl G. Jansky Very Large Array 
(hereafter VLA). The recently enhanced VLA al¬ 
lows us to trace radio continuum emission at levels 
fainter than previously possible via its wide band¬ 
width capabilities. Moreover, observing on a va¬ 
riety of angular scales, the VLA provides distinct 
advantages in understanding disk-halo and halo 
features. For example, both faint diffuse emis¬ 
sion and distinct filamentary structures can be in¬ 
vestigated through the combined use of compact 
and extended configurations. The low frequencies, 
chosen because of their sensitivity to synchrotron 
emission, in combination with observing all polar¬ 
ization products, enables us to derive information 
about the halo magnetic fields and cosmic rays. 


Our thirty-five edge-on galaxies have inclina¬ 
tions higher than 75 degrees. They also adhere 
to limitations on declination (more than —23 de¬ 
grees in order to be observed with sufficient uv 
coverage with the VLA), as well as optical diam¬ 
eter (4'< d 25 < 15') and flux density {Siaghz > 
23 mJy). Three galaxies (NGG 5775, NGC 4565 
and NGC 4244) just outside of these criteria were 
included in the sample as well, due to evidence for 
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extra-planar gas and availability of good ancillary 
data. We refer to Table 1 of Irwin et al. (2012) 


(Paper I) for details of the galaxy sample. 

The full project, with its motivation and sci¬ 


ence goals, is presented in detail in Irwin et al. 


(2012) (Paper I). Two other papers, Irwin et al. 


(2012) (Paper II) and Irwin et al. (2013) (Pa¬ 


per III), present detailed results of CHANG-ES 
observations of NGC 4631 and UGG 10288, re¬ 
spectively. 


In this, the fourth CHANG-ES paper, we 
present all observations which were carried out 
in the shortest baseline array configuration, D, 
and display each galaxy of the survey with its re¬ 
sults in the appendix. In particular, we show the 
Stokes I maps and spectral index maps for each of 
the galaxies in the two frequency bands, as well 
as the polarization map with apparent B-vectors 
superposed, as derived from the Stokes Q and 
U maps. These data products are available for 
download at http://www.queensu.ca/changes 
Additionally, new star formation rates and flux 
densities are presented. 


The paper is organized as follows: in Section 2 
we give a description of the sample selection, the 
setup of the observations and observation details. 
Section 3 describes the data reduction with cali¬ 
bration procedures and Section 4 presents the re¬ 
sulting data products, which are displayed in the 
appendix. Analysis is presented in section 5 and 
the conclusions can be found in Section 6. 


2. Observation setup 

2.1. Survey design 

The details of the observations are presented in 
Table The observations in G-band (central fre¬ 
quency 6.000 GHz) cover a bandwidth of 2 GHz 
(4.979 - 7.021 GHz), in 16 spectral windows and 
1024 spectral channels (64 in each spectral win¬ 
dow). The L-band (central frequency 1.575 GHz) 
observations cover a bandwidth of 512 MHz (1.247 
- 1.503 GHz, 1.647 - 1.903 GHz), in 32 spectral 
windows, and 2048 spectral channels. In L-band, 
we placed the two base bands, of 16 spectral win¬ 
dows each, 144 MHz apart in order to avoid a re¬ 
gion of particularly strong and contaminating ra¬ 
dio frequency interference (RFI). Note that due 
to flagging during the reductions, the final cen¬ 


tral frequencies shown in Tables and for each 
galaxy are adjusted from the values given here. 


2.1.1. Sensitivity 


The theoretical noise level is approximately 6 
/rjy/beam in C-band (with confusion limit and 
flagging taken into account). Due to a high uncer¬ 
tainty in confusion, the theoretical estimate for L- 
band of 89 ^Jy/beam is substantially higher than 
the actual value we achieved, which is around 30- 
35 /rjy/beam. However, due to variable noise in 
the L-band, we measured the noise in areas far 
from the source where it was more uniform (see 
Sec. 4.1). Additionally, a few galaxies exhibit 


higher values for reasons outlined later in this pa¬ 
per. 


2.1.2. Distances 


We have made a few modifications to the dis¬ 
tances used in Table 1 of Paper 1. Five galaxies 
have had their distances modified, whereas the rest 
remain the same as listed in Paper 1. 


The five galaxies with modihed distances are: 
NGC 891, NGC 4244, NGC 4565, NGC 4631, and 
NGC 5907. These were adjusted because they 
have distances derived from the Tip of the Red Gi¬ 
ant Branch (TRGB) |Radburn-Srnith et al.|20TT |. 
Because of the reliability of the TRGB method, 
and that it likely provides the best distances for 
edge-on systems, we have adopted those distances 
here. The most significant distance change was 
for NGC 4565, which went from a Hubble Flow 
distance of 27.1 Mpc to a TRGB distance of 11.9 
Mpc. For the other four galaxies the changes were 
minor (see Table [^. 

Of the remaining galaxies, the distances of the 
five Virgo cluster galaxies (denoted with V in Ta¬ 
ble were adopted from|Solanes et alT (2002). 


All other galaxy distances were taken from 
NED, using the ’’Hubble flow corrected for Virgo 
and the Great Attractor” (HF) distance, and a 
Hubble constant of 73 km/s/Mpc. Since a Hub¬ 
ble flow method may not always be accurate for 
nearby galaxies at low heliocentric velocities, we 
reviewed the galaxy distances for this paper. We 
compared the HF distance with the median of var¬ 
ious Tully-Fisher (TF) derived distances. The lat¬ 
ter are listed in NED as ” red-shift independent” 
distance derivations. Note that we excluded the 
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’’Tully 1988” values listed in NED since those are 
Hubble Flow based. For the majority of the galax¬ 
ies, the changes between the two methods were 
well within the uncertainties associated with ei¬ 
ther method. In particular, it is not uncommon to 
find a factor 2 difference in the various TF-method 
published distances. We therefore retained the 
distances used in Paper I for these 25 systems. 
The comparison showed that for 16 galaxies the 
agreement between the adopted HF distance and 
the median TF distances was better than 25%, for 
9 galaxies the difference was larger than 25%. The 
galaxies with the largest distance uncertainty (> 
35%) include NGC 2683, NGC 3432, NGC 4666, 
NGG 5775, and NGG 5792. See Table 0 for our 
current distance list. 

2.2. Observations 

Of the 405 hours which were awarded for the en¬ 
tire GHANG-ES project, 65 hours were set apart 
for the D-configuration observations in two fre¬ 
quency bands, L and G. The observations were di¬ 
vided up into 13 scheduling blocks, each of which 
contained scans of one primary gain and phase cal¬ 
ibrator (hereafter referred to as the primary cali¬ 
brator) and one zero polarization calibrator to cal¬ 
ibrate polarization leakage from the instrumenta¬ 
tion. Additionally, complex gain calibrations were 
performed using a source (hereafter the secondary 
calibrator) less than 10 degrees from the target 
galaxy. The secondary calibrator was observed be¬ 
fore and after the galaxy scans, every 20-25 min¬ 
utes. The bulk of the data was observed during 
December 2011. Two scheduling blocks were re¬ 
observed in March 2013. See Table [T] for observa¬ 
tion specifics. One scheduling block observed in 
C-band, ID 5062559, suffered increased antenna 
system temperatures at times coinciding with bad 
weather (snow). Gonsequently, the rms noises 
for galaxies NGG 4302, NGC 4388, NGC 4438, 
NGC 4594, NGC 4666 and NGC 4845 increased 
to approximately twice the theoretical expectation 
at C-band. 

2.2.1. Large galaxies 

Eight galaxies in the sample are too large to fit 
inside the primary beam of 7.5' FWHM at C-band. 
We thus observed these in two pointings, sepa¬ 
rated by half the diameter of the primary beam 
at half maximum, i.e. one quarter beam from the 


galaxy centre along the disk on either side of the 
centre. The two-pointing galaxies are marked with 
an asterisk beside their names in Table [T] Table [2] 
lists the coordinates for the separate pointings. 
One single pointing was sufficient for all galaixes in 
L-band, where the diameter of the primary beam 
is 30'. 

3. Data reductions 

3.1. Calibrations 

All calibrations and data reductions were per¬ 
formed using standard routines in the Com¬ 
mon Astronomy Software Applieations package, 
CASAQ using versions 4.2 and earlier. Paper III 
describes in detail the process we have been follow¬ 
ing for calibrations and reductions. A summary 
follows here: 

In order to reduce ringing across the band from 
strong RFI, it was necessary to Hanning smooth 
our data. Antenna-based delay calibrations have 
to be applied prior to this smoothing, which was 
done in an initial stage of the calibrations. 

The absolute flux density scale of the data 
was determined by observing one of two stan¬ 
dard primary calibrators; 3C286 or 3C48 (the lat¬ 
ter for three of our early right ascension galaxies, 
NGC 891, NGC 660 and part of the NGC 2613 ob¬ 
servations), and applying the Perley-Butler 2010 
flux density scale. We used the primary calibra¬ 
tor to determine the bandpass corrections as well. 
See Table [T] for a list of the calibrators which were 
used for each galaxy in the sample. 

3.1.1. Polarization calibration 

For the polarization calibrations, the primary 
calibrator was used as a polarization angle cali¬ 
brator. One of the NRAO standard zero polariza¬ 
tion calibrators was used to calibrate instrumental 
polarization leakage. We corrected for cross-hand 
delays, leakage terms and polarization position an¬ 
gles. Where in some cases the polarization leakage 
calibrator could not be used (due to loss of scan 
during observation or being more heavily flagged 
than the other calibrators), we instead utilized the 
secondary calibrator where the range in parallac¬ 
tic angle was sufficient, i.e. more than 60 degrees. 


^http://casa.nrao.edu 
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This was for example done for the galaxies in a C- 
band scheduling block (ID: 4809749), for which the 
polarization leakage calibrator scan was lost. Q 
and U fluxes were calculated in the standard way 
- see Paper II and references therein. Our tests 
have shown no difference in the result whether the 
secondary calibrator or polarization leakage cali¬ 
brator was used. 

3.2. Flagging 


Additionally, the clean algorithm enabled si¬ 
multaneous htting of a spectral index across the 
band width (’in-band spectral index’) via a sim¬ 
ple power law. In order for this to take place, 
we set the number of Taylor coefficients used to 
model the sky frequency dependence,nterms, to 2 
for all CHANG-ES galaxies. We defer the details 
of spectral information to Section [3^ and will de¬ 
scribe the spatial fitting, as implemented for the 
CHANG-ES project, in this section. 


For all scans, calibrators and galaxy scans alike, 
we inspected the uv-data by eye. All bad data, 
whether from RFI or instrumental effects, were 
manually flagged. A number of tests were car¬ 
ried out using automatic CASA flagging routines 
as well as the CASA pipeline; however, we found 
significantly improved results by continuing with 
manual flagging. 

Generally, fewer antennas could be used for our 
L-band observations than for C-band, since not all 
VLA antennas had yet been equipped with L-band 
receivers at the time of observation. Additionally, 
there is typically more RFI at L-band. 

Calibrations and flagging were carried out iter¬ 
atively, until a well-calibrated data set could be 
obtained. The final calibrated galaxy scans were 
then separated from the calibrator scans, in order 
to be imaged in Stokes I, Q and U. 


3.3. Imaging 


We have produced Stokes I, Q and U images, 
as well as spectral index maps (a, where /jy oc ^'“) 
and uncertainty maps, Aa (see Sect. 3.41 from the 
calibrated data. 

The clean algorithm in CASA, which we use 
for our wide-held synthesis deconvolution, allows 
for multi-scale multi-frequency synthesis, ms-mfs 
(see Rau & Cornwell 2011, for a full descrip¬ 


tion). We set the multi-scale feature to look for 
hux components at a variety of spatial scales, to 
account for the extended emission of our galax¬ 
ies. Typically our scales ranged from zero (’classic 
clean’) up to approximately hve times the synthe¬ 
sized beam, but sometimes required adjustments. 
In L-band, and for the uv-tapered weighting in 
particular, fewer scales were at times necessary 
in order to remove artifacts for point-like sources; 
for example a classic clean was set for NGC 2683, 
NGC 2992, NGG 4157 and NGG 4845. 


3.3.1. Weightings 


For each galaxy, the Briggs robust weighting 


(Briggs 19951 was used for both G and L band 


data. Additionally, we imaged a uv-tapered ver¬ 
sion of this weighting, which we introduced in or¬ 
der to emphasize the broadscale structure^ We 
achieved an increase in beam size of approximately 
80% in C-band and 40% in L-band, by typically 
applying a 6 kA Gaussian taper in C-band and 
a 2.5 kA Gaussian taper in L-band to the uv- 
distribution. Image products of both these weight¬ 
ings are included in Data Release 1 and are re¬ 
ferred to as robust 0 and uv-tapered weightings, 
respectively. 


In G-band, images sufficiently clear of artifacts 
with uv-tapered weighting could not be achieved 
for NGG 4438 and NGG 4845 due to particu¬ 
larly strong contaminating sources in the field. 
In L-band, the uv-tapered weighting was either 
challenging to achieve (for the same reason) or 
unnecessary (the sources were imaged as point 
sources already at robust 0) for six of our galax¬ 
ies; NGG 3448, NGC 4192, NGC 4388, NGC 4438, 
NGC 4845 and UGC 10288, and the attempt was 
discarded for these. Hereafter any images shown 
will be of robust 0 weighting, unless otherwise 
specified. 


3.3.2. Stokes I imaging 

Each imaging run was cleaned down to a flux 
density level of 2.5-3 a. With few exceptions, 
we imaged the entire field (i.e. without specify¬ 
ing regions), since a lower rms could be obtained 
with this strategy. Additionally, the vast major¬ 
ity of our data do not suffer from significant arte¬ 
fact patterns with spurious sources outside the 


^For an example of a different weighting, see Paper II for 
which a naturally weighted map of NGC 4631 was included. 
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source region, which could jeopardize this strat¬ 
egy. When needed, one to two self calibrations 
were performed in order to deal with remaining 
artifacts. 


Our self calibration strategy was based on iter¬ 
atively refining the calibration table, by self cali¬ 
brating at successively deeper thresholds and thus 
improving the input model, which in turn acts on 
the original visibilities. At each step, care was 
taken to a) check the model to be used for the self 
calibration that it did not include any artifacts 
and b) check that the peak intensity of the source 
would not decline (such a scenario was aided by 
opting for a phase-only self calibration at the first 
iteration, or by applying the self calibration at a 
shallower threshold). 


In some cases, our cleaning strategy had to 
be adapted to the circumstances of the data (the 
quality of a dataset, strength of emission, contam¬ 
ination in the field, etc). This particularly ap¬ 
plies to data with high dynamic ranges (HDR), 
due to galaxies with strong, often point-like cen¬ 
tre sources or strong field sources (for example 
NGC 660, NGC 4594, Virgo galaxies). Our strate¬ 
gies to deal with these situations include one or 
more consecutive runs of self calibrations, speci¬ 
fying regions, attempts at peeling for strong field 
galaxies (see |3.3.6 ) and clean parameter adjust¬ 
ments (including multi-scale adjustments). Gen¬ 
erally, the best one can achieve in terms of HDR, 
is a signal over noise of the order of 10000:1 (S. 
Bhatnagar, priv.comm.). Our dynamic ranges are 
listed in column 8 of Tables H] and [5l Addition¬ 
ally, in cases of strong centre sources, choosing a 
smaller cell size will characterize the steep beam 
better and render a cleaner result (e.g. NGC 4594 
and NGC 4845, both in C-band). 


3.3.3. Polarization imaging 

Stokes Q and U images were produced similarly 
to I. These maps were used to create polarization 
intensity and polarization angle maps. The po¬ 
larization intensity maps were corrected for bias 
caused by the noise levej^and the polarization an¬ 
gle maps are cut off at the Scr level. We note that 


®Note that the maps shown in panels (e) and (f) of the 
images in the appendix are the noise biased maps. Peak 
polarization values have however been measured from the 
corrected maps. 


if the fraction of polarization intensity over Stokes 
I intensity is less than 0.5%, the polarization is not 
believable since the calibration cannot be guaran¬ 
teed below that level (S. Myers, priv. comm.). 

A uniform Faraday screen resulting from the 
ionosphere is corrected for in the standard polar¬ 
ization data reduction. However, if there is differ¬ 
ential Faraday rotation, for example between the 
location of the primary calibrator and the source, 
then this is not corrected foi[^ This effect is negli¬ 
gible at C-band. There could, however, be minor 
effects at L-band, which will be investigated in a 
subsequent polarization paper. 

Our vector maps (panels (d), (e) and (f) of 
the figures in the appendix) are also ’apparent B- 
vectors’ which simply are E-vectors rotated by 90 
degrees, because internal Faraday rotation in the 
galaxy has not been corrected for. 


3.3.4. Primary beam corrections 


Wideband primary beam corrections were car¬ 
ried out with the CASA task widebandpbcor af¬ 
ter cleaning. This task accesses information about 
the beam from the visibility data, and calculates 
a known NRAO-supplied model of the primary 
beam (PB) at each frequency channel in each 
band. A PB cube is formed, and Taylor terms are 
found in the same way as was done for the data 
set when imaging. Each Taylor coefficient image 
(two in our case, see Sect. 3.3 and 3.4.1) is then 
corrected. 


This task also corrects the spectral index im¬ 
age above a given threshold (we use 5cr), as well 
as corrects a map of the formal error of the spec¬ 
tral index (see Sec. 3.4.2). Subsequently, we also 
applied primary beam corrections to the polariza¬ 
tion intensity maps, by dividing the image with 
the beam created in the widebandpbcor task. 


Note that in this paper, any displayed images 
are not PB corrected, but any measurements are 
made from the PB-corrected versions of the maps, 
unless otherwise indicated. For total intensity im¬ 
ages, use of a model for the PB rather than an 
observation-specific PB introduces negligible error 
over the scales of our galaxies. However, the effect 
on the spectral index maps must be investigated 
in more detail, which we do in Sec. |3.4.2[ 


^The software is not yet available in CASA. 
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3.3.5. Mosaicking two pointings 

For the eight large galaxies which were observed 
in two pointings at C-band (see Table , each in¬ 
dividual pointing was calibrated and imaged sep¬ 
arately. Although it would be preferrable to mo¬ 
saic the pointings onto the same grid while still in 
the uv-plane, this feature is not yet implemented 
in CASA at the time, for the clean settings that 
we required (i.e. nterms = 2). Instead, we com¬ 
bined the images from the two pointings, forming 
a weighted (by the primary beam) average in the 
overlapping region. The images were combined 
in the overlapping regions out to 10% of the pri¬ 
mary beam. This rather inclusive choice of cut-off 
worked well for the D configuration data, but will 
likely be adjusted to 50% for the B and C config¬ 
uration data in the future. 

This strategy was applied to both the non-PB- 
corrected images, as well as the two primary beam 
pointings. The combined Stokes I image was di¬ 
vided by the combined primary beam in order to 
do the PB correction. 

The two-pointing spectral index maps, which 
had already been individually PB-corrected, were 
similarly combined, as well as their accompanying 
error maps. 

We assessed whether the two pointings had dif¬ 
ferences within the overlap regions, such as in 
noise levels. In such cases additional weightings, 
based on the noise level, would have to be consid¬ 
ered. This was however not the case for any of 
the CHANG-ES galaxies, where separate point¬ 
ings generally had been observed in turns in the 
same scheduling blocks, or in similar conditions. 


3.3.6. Peeling of strong field sources 


For a couple of cases, for example NGC 660 and 
NGG 4388, bright field sources interfere strongly 
in the cleaning process, causing a less than ideal 
imaging result fraught with severe artifacts (L- 
band only). An attempt was made at peeling these 
interfering sources in those examples, that is, to 
remove or decrease the intensity of the interfering 
source so that its impact is lessened in the final de- 
convolution. The process includes self calibrations 
with the interfering source centred, and subtract¬ 
ing the source from the uv-data using the GASA 
task uvsub, in several steps. For more informa¬ 


tion on peeling, see Pizzo & de Bruyn (2009), as 


well as the PhD thesis by Adebahr (2013). Table 


lists the sources we peeled from NGG 660 and 
NGG 4438. 


3.4. In-band Spectral Indices 
3 . 4 . 1 . Formation of Spectral Index Maps 


As indicated in Sect. |3.3| the ms-mfs algorithm 
allows for spectral fitting across the wide bands 
that have been used in the GHANG-ES survey. 
Thus, a single observation over a single wave-band 
permits the determination of spectral indices, a, 
as a function of position within the galaxy. As 
described in Rau & Gornwell ( 2011[ ) (see also Pa¬ 
per II), a specihc intensity, 1^, at some frequency, 
V, within the band, can be fit with a function of 
the form ^ where /3 is a cur¬ 

vature term and vq is a reference frequency (our 
central frequency, vq, given in Tables and [^. As 
implemented, such a htting function is expanded 
in a Taylor series about t'oj such that the first Tay¬ 
lor term (the map, TTO) corresponds to a map 
of specific intensities at vq and the second Tay¬ 
lor term (the map TTl) contains information on 
a such that a = TTl/TTO. The third Taylor 
term (TT2) allows for the recovery of the spectral 
curvature, (3, where /3 = TT2/TT0 — a{a — l)/2. 

In principle, any number of Taylor terms could 
enter into such an expansion with increasing num¬ 


bers of terms improving the spectral fit (see Rau 


& Gornwell 2011, for examples). In practice, how¬ 


ever, the effective number of terms is limited by 
the signal-to-noise (S/N) ratio and calibration pre¬ 
cision of the data. For the GHANG-ES galaxies, 
tests that we have run gave poorer results with 
3 terms than with 2; that is, the rms noise in the 
images is higher for such fits, when cleaning to the 
same threshold. 

The curvature maps also show large and some¬ 
times unrealistic variations, point-to-point. While 
globally averaged values of /3 can still be useful 
(see Paper II for an example), we have adopted 
two Taylor terms for the GHANG-ES project as 
indicated earlier. A two-term Taylor fit then takes 
the simple form {{v — Vff)/Vff). 

A map of random errors describing the accuracy 
of the fit, A a is also formec0 

The beauty of the spectral fitting is that it is 


®See Eqn. 39 of Rau & Cornwell (2011 1 . 
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carried out by flux component which means that 
a uniform resolution is achieved over the band, 
with a restoring beam applied to the TTO and TTl 
maps at the end of the fitting process. Also, the 
S/N of the entire band applies to each flux com¬ 
ponent. If one were to form a spectral index in 
the classical way within a band, one would need 
to break up the band into its various channels (or 
groups of channels), smooth each channel to the 
poorest resolution of the lowest band frequency, 
and accept the very high noise of these channels 
as a cut-off when forming a. maps. 

Clearly, the method, as implemented in CASA, 
is far superior for any individual banc^ Neverthe¬ 
less, there are limitations, as we outline below. 

3.4-2. Post-Imaging Corrections to the Spectral 
Indices 

Because the primary beam (PB) varies with fre¬ 
quency, it imposes its own spectral index, apB, 
onto the a maps. This effect is significant; for ex¬ 
ample, at the 50% level at L-band, apB = —1.6. 
We use the CASA task widebandpbcor, described 
in Section [3.3.4[ to correct for aps, and at this 
point, we cut off the a maps at 5a, where a is the 
rms of the image map (i.e. the TTO, or map). 

Once PB-corrected a maps are formed, they 
can show large variations around their peripheries 
and the corresponding error maps, A a, show cor¬ 
respondingly large errors. In addition to the 5cr 
cut originally applied, we also cut off all a maps 
wherever the corresponding error map exceeds a 
value of 1.0. The choice of this cutoff is arbitrary; 
however, various tests showed that more aggres¬ 
sive cuts (e.g. where A a > 0.75) resulted in 
undesirable ‘holes’ in some a maps. Note that 
this step does not change the PB correction of the 
a maps. It simply discards data points that are 
known to have large errors. 

Both a and A a maps are computed for each 
map pixel (cell). However, these pixels are not in¬ 
dependent. For example, at L-band, there may 
typically be 64 cells/beam and at C-band, 50 
cells/beam (both uv-tapered). Consequently, the 
‘per pixel’ errors are much greater than the ‘per 
beam’ errors and the a maps, as produced, show 
variations pixel-to-pixel which are larger than 


®For an example showing a result when more than one band 
is used, see Paper III. 


beam-averaged variations would be. In addition, 
maps of A a show artifacts, i.e. regions of low 
A a that have a ‘thread-like’ appearance through 
the map. These artifacts depend on the choice 
of multi-scales that are used during imaging and 
shift position if the map is remade with different 
scales. If the adopted scales are changed, these 
’threads’ also shift position. 

To ameliorate the above issues, we have ‘smoothed/ 
averaged’ our maps of a and A a over the size of 
a beam. This is accomplished by artificially in¬ 
serting into the map header a Gaussian beam size 
whose area (Gaussian-weighted) is equivalent to 
the pixel area and then convolving the map to 
the correct clean beam siz^H The result does not 
change the spatial resolution of the images but 
recognizes the non-independence of the pixels and 
strongly minimizes pixel-to-pixel variations. For 
example, the rms of an a map could decrease by 
as much as a factor of 2 after such a convolution. 
These are the final maps that are shown in the 
appendix, panels h), i), k) and 1). 

In summary, our final spectral index maps and 
all calculations performed on them apply to a after 

1. PB-correction, 

2. a 5(7 cutoff, 

3. a cutoff wherever A a > I.O, and 

4. convolving with a Gaussian to smooth over 
pixel-to-pixel variations within any clean 
beam. 

3 . 4 . 3 . Uncertainties in Spectral Index Maps 

As indicated above, a is determined via fitting 
flux components over the band during the imag¬ 
ing process. Fitting errors should, for the most 
part, be accounted for by the A a maps, but those 
maps consider random errors only. Therefore, a 
variety of tests have been carried out on the in- 
band spectral index maps to help us understand 
realistic uncertainties, and we outline them here. 
Largely qualitative descriptions are provided, with 
a quantitative summary at the end. In the fol¬ 
lowing, we refer to the final a maps, as described 
above, unless otherwise indicated. 

Edge-effects: 


^The CASA task imsmooth is used. 






In spite of the cutoffs and smoothing outlined 
in the previous section, values in the spectral in¬ 
dex maps (as well as their errors) tend to become 
extreme at the edges (for example, see the L-band 
panels (i) and (1) of NGC 3735 in the appendix). 
Primary Beam Correction: 


The NRAO-supplied PB model cannot take 
into account variations in the PB that are spe¬ 
cific to any given observation. The true PB will 
change with flexure in the antennas and therefore 
with position on the sky. It will also rotate on the 
sky as the source is tracked and therefore any non- 
axisymmetry in the PB will also introduce erroij^ 
As a result, PB correction errors, which increase 
with distance from the pointing centre, must be 


understood. Bhatnagar et al. (20131 have shown 


that errors in the a maps are negligible out to the 
half-power point but increase significantly beyond 
that distance (see their Fig. 4, right). For L-band, 
our largest a map (NGC 4565) extends to a radius 
of 6.8 arcmin which is well within the HPBW of 
the L-band PB (radius of 15 arcmin, or FWHM of 
30 arcmin). At C-band, however, the PB radius 
is only 3.75 arcmin and consequently a maps with 
radii larger than this will show such errors in their 
outer regions. 


We have investigated such potential PB errors 
in two ways: 

The first is to examine a maps using cases in 
which the observations of a given galaxy were car¬ 
ried out in two different observing sessions, such 
as NGC 660 (see Table [^. As such, it is possible 
to form two different a maps for the same galaxy 
and same sky pointing but observed on different 
days. This has allowed us to compare the result¬ 
ing a maps and to compare the variations between 
those days with the random errors as given in the 
A a maps. We tested both low and high dynamic 
range cases and compared the two results quanti¬ 
tatively. 


The second is to examine our C-band observa¬ 
tions of large galaxies where two offset pointings 
were observed (see Sec. 2.2.11. A comparison of 
a made for the individual pointings allowed us to 


^Recently, it has become possible to correct for narrow- 
band time-varying PB effects during the imaging process 
using the ‘A-Projection’ algorithm (|Bhatnagar et al. |2008 
|2013[ l. However, an algorithm applicable to CHANG-ES 
data (i.e. for wide bands) is not yet available or practical. 


investigate a errors at larger distances from the 
PB centre. 


As an example, Fig. (panel a) shows a map of 
the absolute pixel-by-pixel difference between the 
PB-corrected spectral index measurements in the 
two C-band pointings of NGC 891. An increase of 
the differences with distance from the map centre 
is clearly apparent. We formed the spectral index 
error (shown in panel b) by computing the PB- 
weighted deviation from the average spectral index 
of the two pointing^ 


This error map is for the most part comparable 
to the error calculated by the ms-mfs algorithm 
(panel c)), as the map of the ratio of these two 
errors (panel d)) illustrates. A major exception 
to this is the disk of the galaxy, but only in those 
parts that are located outside the 70% PB level 
(i.e. where the PB gain is less than 0.7) of either 
pointing. The average error ratio in these narrow 
regions around the mid-plane is ~5 (with individ¬ 
ual pixel values up to ~15), i.e. here the error orig¬ 
inating from the pointing differences is on average 
~400% higher than the error resulting from the 
ms-mfs spectral fitting (for comparison, the aver¬ 
age error ratio of NGC 4565 is 3.3). Such high 
ratios are primarily a consequence of the small 
ms-mfs-based errors in regions of high signal-to- 
noise, yet the increase of the pointing differences 
with distance from the pointing centres is signif¬ 
icant. In particular, beyond the half-power point 
the mid-plane errors in panel b) of Fig. [^increase 
up to ~0.5, whereas in panel c) they remain be¬ 
low 0.1 throughout the disk, as shown by the dis¬ 
played contours. While the errors in panel b) are 
in rough agreement with iBhatnagar et al. (2013) 


in the sense that they do not exceed 0.1 out to ap¬ 
proximately the half-power point (not considering 
the above-mentioned edge effects), the error ratio 
map suggests that the simple PB model used by 
the widebandpbcor task already shows significant 
inaccuracies at the 70% level. 


Resolved out structures: 


Although the method of calculating the spec¬ 
tral index ensures that the resolution is common 
across the band, if there are structures that are 
resolved out at one end of the band compared to 
the other, this is not accounted for in the spectral 


® - Oavg)'^ + W 2 {a 2 - otavg)'^)/+W 2 ), where 

wi = beamilbeam 2 and W 2 = beara 2 lbeam\ 
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index maps. For example, there may be struc¬ 
tures resolved out at the high frequency end of 
the band, but not at the low frequency end, which 
would steepen the spectral index. 


Comparison with Classical Spectral Index: 

A comparison was made between the in-band 
spectral index and a classically formed spectral in¬ 
dex from a given band. The classically formed a 
maps were made from the spectral index maps af¬ 
ter the PB-correction was applied but not after ap¬ 
plying further processing (i.e. PB-correction and 
a 5a cutoff were applied, but not a cutoff based 
on the error map, nor smoothing to the size of the 
beam, as described in Section 3.4.2); they were 
formed by imaging the lowest end of the band 
and the highest end of the band, smoothing to 
a common spatial resolution and then combining 
in the classical way. The results were in agreement 
within errors. 


Two-pointing PB cutojfs: For ‘two-pointing’ 
galaxies (large galaxies at C-band) the maps, as 
indicated in Section |3.3.5[ were combined over 
regions in which the value of the primary beam 
exceeded 0.1 (10% of the peak). When mosaicking 
is carried out, each point is weighted by the PB 
such that points that are farther from the beam 
centre are weighted down and therefore the differ¬ 
ence in total intensity maps, whether one cuts off 
points where the PB is > 10% or where the PB 
is > 50%, is entirely negligible. However, since a 
maps are known to increase quite strongly below 
the 50% level, we remade our mosaicked a maps 
with a 50% cutoff to compare it to the 10% cut-off 
values used in this paper. For these results, we 
found only minor differences (for example, if the 
two maps are subtracted, the rms in the regions 
between the 50% and 10% PB cutoffs is < 0.1). 


d) The largest source of uncertainty in A a relates 
to the PB model correction away from the point¬ 
ing centre, which most strongly affects our largest 
galaxies in C-band. Spectral index errors increase 
significantly with distance from the pointing cen¬ 
tre, such that beyond the 70% PB level of either 
pointing (i.e. at distances greater than 2'25" from 
the respective pointing centre at C-band) the A a 
maps typically underestimate the true errors in 
the disk by a factor ^5 (see Fig[^. 


4. Results and data products 


A wide range of data products are created for 
each galaxy. We have assembled a selection of 
these for each galaxy, which are presented in the 
appendix. See the appendix for a detailed listing 
of each panel. In summary, the 12 panels show: 
(a) and (b): the Stokes I maps for 2 C-band 
weightings (robust 0 and a uv-tapered weight¬ 
ings), 

(c) : Stokes I map for one L-band weighting (ro¬ 
bust 0 only), 

(d) : an optical image with contours and apparent 
B-vectors, both from the uv-tapered weighting of 
C-band 

(e) and (f): polarization intensity maps with ap¬ 
parent B-vectors overlaid for C and L-bands, re¬ 
spectively 

(g) : a composite image of different weight¬ 

ings/bands (see 4.2), 

(h) , (i), (k) and (1) spectral index maps (h, i) with 
corresponding error maps (k, 1) in C and L-bands 
respectively. 

(j) a wide view of the L-band field. 


Tables and list beam sizes, rms noises and 
dynamic ranges for both bands and weightings. 


In summary, our final a maps should present 
realistic results with the following cautions: 

a) Extreme values around the edges are artifacts 
and should be ignored. 

b) Measurements of a should not be quoted ‘per 
pixel’ but rather averaged over a beam when car¬ 
rying out evaluations or comparisons. 

c) The calculated A a maps typically underesti¬ 
mate the true errors by « 20 % for galaxies of small 
angular siz^^ 

^'^Estimated from results where galaxies have been observed 


4.1. Rms noise 

All rms noise values mentioned throughout 
this paper have been measured from non-primary 
beam corrected maps. For C-band, the rms noise 
is measured as an average of the regions through¬ 
out the imaged field of view which do not contain 
detectable background sources. 

For L-band, however, there are so many back¬ 
ground sources near the galaxy that it is difficult 


in two separate observations. 
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to find sufficiently large regions within which the 
rms can be consistently measured. Residual clean¬ 
ing artifacts also tend to be larger at L-band, es¬ 
pecially close to the galaxy where the primary 
beam response is high. Therefore the rms tends 
to be variable when measured near the source, 
but declines to more consistent values with dis¬ 
tance from galaxy as the primary beam response 
declines. For consistency, we quote rms values far 
from the galaxies where there are smaller varia¬ 
tions in rms between measurement regions. Rms 
levels near the galaxy could be up to a factor of 
two higher than the quoted values of Tables and 
but it is apparent from the panel c) of the ap¬ 
pendix figures, that our choice of 3 ct as the lowest 
contour generally well represents the faintest be¬ 
lievable emission. 

4.2. Displaying the halo 

Panel (g) of the figures in the appendix pro¬ 
vides an exploratory image intended to help the 
viewer discover information about individual ra¬ 
dio halos. For example, one may like to know if 
point sources in the disk (anti-)correlate with very 
extended halo emission. The CHANG-ES survey 
is data rich, with two (or more) weightings pro¬ 
vided per band and overlaying contour plots of 
these for both bands can result in confusing dia¬ 
grams. Therefore we overlay colourized ’’transpar¬ 
ent” images of the weightings so that the viewer 
may relatively quickly apprehend the relationship 
between structures in the two observing bands. 
Additionally, artifacts, confusion structures, and 
non-random noise are difficult to remove mathe¬ 
matically from radio data. Our visualization ap¬ 
proach uses masking in order to mitigate ’’back¬ 
ground” artifacts. 

First the fits data are stretched using the 
KARMAp^ visualization package’s kvis task. On 
intensity-inverted, logarithmically scaled data, we 
use the ’’greyscale 3” option in the pseudo-colour 
tool for adjustments. The resulting images are 
saved in eps format and used as input into the Gnu 
Image Manipulation Program (GIMPj^ This 
package allows the user to stack images in ” layers” 
that can be combined, as if they are transparent, 
using a variety of blending mode algorithms. Usu- 


^^See http://www.atnf.csiro.au/computing/software/karma/ 
^^See gimp.org 


ally we stack four images, i.e. two weightings of 
each observed band. We colourize the G band data 
(blue for the robust 0 weighted data and green for 
the uv-tapered weighting) while the L band data 
are left as greyscale images. The order in which 
the images are combined, and which algorithms 
are used, are described in Fig. In the resul¬ 
tant qualitative image, often both point sources 
in the disk and diffuse halo structure are evident 
simultaneously, s 

4.3. Comments on individual galaxies 

A few galaxies warrant some extra comments, 
either because of differences in the data reduction 
procedures, or because the results were interesting 
or unusual and therefore caught our attention. In 
this section, we list these galaxies with comments. 
Note that this is not meant to be a thorough dis¬ 
cussion of each galaxy. 


4.3.1. NGC 660 


The data of polar ring galaxy NGG 660 have 
the highest dynamic range in both bands, due to 
a strong central source. Gonsequently, the result¬ 
ing images have higher than expected rms values, 
particularly in L-band. Nevertheless, the achieved 
signal to noise of the G-band map with robust 
weighting reached 70000, on the higher end of 
what best can be expected (see Section [3. 3. 2[ ). 

In spite of careful cleaning, including peelinj^ 
of the nearest of the two strongly interfering field 


sources in L-band (see Section 3.3.6), artifacts are 


still present. This could potentially be affecting 
the spectral index results, where regions of flat 
spectral index are crossing the disk. However, 
these regions could also be an effect of the pos¬ 


sible AGN (see 5.2.1). 

We note that the polarization intensity is low, 
with its fraction of polarization over Stokes I in¬ 
tensity just a 10th of the 0.5% which is considered 


a believable signal (see Section 3.3.3). 


4.3.2. NGC 3556 and NGC 5775 

Both of these galaxies are well-known to have 
significant halos, which we can also see in our data. 


^^Only the Stokes I and spectral index images were produced 
from the peeled data, while polarization images were not, 
since the off-centre confusing sources are not as problematic 
in Q and U as for total I. 
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Moreover, they both show very flat spectral in¬ 
dices at C-band throughout the disk, where the 
errors in spectral indices are the lowest. Neither 
data set showed any particular problems and were 
single pointings. Bearing in mind the summary 
regarding spectral index in Section 3.4.3[ this sug¬ 
gests that thermal emission may be more strongly 
dominant for these galaxies in their disks. 

For NGC 5775 at L-band, some broadscale po¬ 
larization features have been found in the field that 
are most likely from foreground Galactic emis¬ 
sion (this is more obvious when a larger field than 
shown in the appendix is displayed). The Galactic 
coordinates of this source place it almost directly 
over the Galactic centre, although at reasonably 
high Galactic latitude. Indeed, the Galactic coor¬ 
dinates of NGC 5775 {I = 359.4°, 6 = 52.4°) in¬ 
dicate that we view it through the northernmost 
tip of the Fermi bubbles that extend 55 degrees 


above and below the Galactic centre (Ackermann 


et al. 20141. Substantial polarized emission at 


2.3 and 23 GHz has been found to coincide with 
the Fermi bubbles, including ridge-like filamentary 
structures crossing through this particular loca¬ 
tion (Carretti et al. 2013). Our polarization im¬ 
ages are likely picking up this same extended fore¬ 
ground structure at lower frequency. 


4.3.3. Three Virgo galaxies: NGC4192, NGC4388 
and NGC 4438 

These galaxies are highly affected by contam¬ 
ination from strong sources in the field, such as 
M87, rendering higher than normal rms noise and 
artifacts hard to clean out. This may have an 
effect on both spectral index results and polariza¬ 
tion. 

We note that the C-band spectral index is steep 
in the inner disk for NGC 4192 (a rather weak 
source at the low end of our flux density cutoff for 
the survey), and that the uncertainties are high. 

NGC 4388 was particularly affected by M87, 
and we made an attempt at peeling the source in 
L-band, as well as self calibrating on M84 (the 
other disturbing source in the field), which ren¬ 
dered some improvement. The resulting rms at 
L-band was too high for us to detect significant 
polarized flux (panel f). 

Also NGC 4438 is strongly affected by residual 
side lobes from M87. Despite attempts of self cali¬ 


bration and peeling, the rms could not be brought 
down to lower than almost three times the ex¬ 
pected value. L-band observations of NGC 4438 
differ from the other galaxies, in that only the up¬ 
per half of the frequency range (spectral windows 
16-31) was used for imaging, with a central fre¬ 
quency of 1.77 GHz. Q and U images were made 
over the same upper half of the frequency range in 
order to match the total intensity image, as well 
as the spectral index image. 


4 . 3 . 4 . NGC 4594 


The strong centre of NGC 4594 (the Sombrero 
galaxy) resulted in residual side lobes which made 
it difficult to detect the weak east-west disk for the 
two-pointing data set in C-band. Careful selection 
of self calibration inputs eventually revealed the 
disk, and subsequent merging of the two pointings 
further helped to cancel out artifacts from the cen¬ 
tre. 


In L-band, a plume is seen to the north of the 
core and is roughly in the direction of the jet ob¬ 
served at much higher resolution by |Hada et al.| 
(2013). An unusual polarization structure and a 
flat spectral index associated with the core and 
towards the north are also observed. A more con¬ 
servative spectral index cut-off could be beneficial 
for this data set (for example, a lOcr cut-off was 
adopted for the in-depth study of NGC 4845, see 
Irwin et al. 2015 (accepted by ApJ). 


4 . 3 . 5 . NGC 4631 

The C-band spectral index shows an unex¬ 
pected flattening from the central to the outer disk 
for this two-pointing data set. This gradual devi¬ 
ation has not been detected in previous observa¬ 
tions of this galaxy (Hummel & Dettmar 1990| ), 
but is within the errors pointed out in point d) of 

Sec. [1331 


4 . 3 . 6 . NGC 4666 

Despite its small angular size well within the 
symmetric primary beam, the C-band spectral in¬ 
dex displays an asymmetry between the two halves 
of the disk, such that it is flatter on the south-west 
side compared to the north-east. This asymmetry 
corresponds to an asymmetry in the polarization 
at the same frequency. Interaction with a nearby 
companion, NGC 4668 (to the south-east) could 
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be a factor, and/or motion through an intergalac- 
tic medium. No other technical issues (antenna 
flexure, solar influence) are problematic. 

13.1. NGC 4845 

It is worthwhile to note that this galaxy dis¬ 
plays an interesting variation in flux density 
level and further exploration, being published in 
CHANG-ES V (Irwin et ah, 2015, accepted by 
ApJ), indicates with little doubt that the source 
is indeed variable. CHANG-ES V also explores 
the spectral index and polarization results in more 
detail. 

4.3.8. NGC 5084 

Due to the low declination of this galaxy, the 
C-band observations were divided up into three 
scheduling blocks in order to a) avoid shadowing 
and b) increase uv-coverage. The scheduling block 
in which the first pointing was observed did not 
have a polarization leakage calibrator scan, and in¬ 
stead the secondary calibrator used for NGC 4192 
was used to calibrate polarization leakage (the sec¬ 
ondary calibrator of NGC 5084 did not have a 
sufficiently large parallactic angle span). Part of 
the second pointing was strongly affected by shad¬ 
owing and only 21 antennas were used for these 
scans. 

The C-band images show emission extending 
east and west of the point-like centre, as well as 
southward from the eastern extension - these are 
likely real, since their intensity is greater than lOcr 
(and they roughly follow the disk). However the 
extensions seen above and below the centre may be 
spurious since they seem to align with weak resid¬ 
ual side lobes. In L-band, we find no apparent 
polarization from the source. The spectral index 
is flat in L-band and flat to slightly positive in C- 
band (but the positive trend coincides with higher 
error values). 

4.3.9. UGC 10288 

UGC 10288 fell within the intensity cutoff of 
the CHANG-ES sample by piggybacking on the 
previously unresolved strong background source to 
the west of the disk. We refer to Paper HI for 
full details on the observations and analysis of this 
galaxy. 


5. Analysis 

5.1. Flux densities 

The flux densities listed in Golumns 2 and 4 
of Table [7] are the average of measurements taken 
at the two weightings (robust 0 and the robust 0 
with a uv-tapering applied). Our measurements 
avoided disconnected point-like features immedi¬ 
ately around the source but did include anything 
that might be considered an ’extension’ (at this 
point it is unknown whether the extension is in¬ 
deed related, or not). For example, NGC 3448 
shows clearly a point source at C-band, but the 
same is not visible as a distinct source at L-band. 
The error (Columns 3 and 6), which is the dif¬ 
ference between the two weightings, also encom¬ 
passes inclusion vs non-inclusion of such point 
sources. In most cases, however, the calibration 
error (flux density scale accuracy) which is esti¬ 
mated to be about 29cp^is larger than this error. 
The error columns in the table cite the larger of 
the two uncertainty estimates. 

For some cases, such as NGC 5907, a back¬ 
ground source is visible through the disk - in such 
cases, the flux density from these sources has been 
removed from the total flux density measurement. 

Our new flux density values are the best so far 
at these two frequencies (cf. Paper I). 


5.2. Star formation rates 


Star formation rates and "surface densities” 


(see below) are based largely on WISE (Wright 
et al. 2010| 22-p,m images with resolution en¬ 
hanced over the standard survey products to a 
final value of 12.4" via the WERGA (WISE En¬ 
hanced Galaxy Resolution Atlas) process (Jar- 
rett et al. 2012 20131. WERGA images in all 


four WISE bands were kindly provided by T. 
Jarrett. Foreground stars were removed through 
PSF subtraction in the WERGA pipeline. Back¬ 
grounds/foregrounds are relatively uniform, al¬ 
though one galaxy showed a significant jump in 
background level at 22 ^m close to the galaxy but 
not crossing it. A constant background was there¬ 
fore subtracted. Flux densities were measured 
with apertures large enough to include all emis¬ 
sion from each galaxy and converted to Jy using 


^^https: / / science.nrao.edu/facilities/vla/docs / manuals / 
oss2014a/performance / fdscale 
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conversions given in the WISE Explanatory Sup¬ 
plement (Cutri et al. 20111. Flux densities from 
any background sources were subtracted. Three 
corrections to the flux densities, as described by 


Jarrett et al. (20131 are applied; an aperture cor¬ 


rection for extended sources, a small colour correc¬ 
tion appropriate for dusty, star-forming galaxies, 
and a calibration correction at 22 /rm appropri¬ 
ate for spiral galaxies. Formal uncertainties on 
the flux densities are small, and the final uncer¬ 
tainties include a 1.5% flux calibration uncertainty 
(Jarrett et al. 20131, although by far the biggest 
source of uncertainty in the SFR is the galaxy dis¬ 
tance. No corrections were made for extinction. 

Spitzer MIPS 24-/xm images ( Rieke et al.|2004 |, 
at a resolution of 5.9" and with a range of sensi¬ 
tivities, were found in the NASA/IPAC Infrared 
Science Archive for 16 of our galaxies, and will 
mostly be used in future work, but are employed 
for one purpose below - hence, their reduction 
is also described here. The ”pbcd” mosaics pro¬ 
duced by the standard MIPS pipeline reduction 
were found to have low-level ’’jailbars” in most 
cases. We therefore reformed the mosaics from the 
individual’’bed” images via ”self-calibration” (di¬ 
viding each bed image by the median of all images) 
with the MOPEX software. This procedure re¬ 
moved the jailbars very well. No other image arti¬ 
facts were noted. Multiple resulting mosaics were 
added to produce final images. Background gradi¬ 
ents were then subtracted with the IRAF program 
background. 

Jarrett et al. (2013) have shown that WISE 


22-/xm fluxes {vF^,) correlate extremely well with 
Spitzer 24-^m fluxes for their sample of 17 nearby 
galaxies, allowing star formation rate calibrations 
developed for the latter by Rieke et al. (2009) to 
be applied to the former. We therefore use the 
relation between 22-/im vLi, and SFR given by 


Jarrett et al. (2013) to calculate SFRs, except for 


NGC 4388 where the Rieke et al. (2009) relation 
is used. 

Fig-i shows the same 22 vs 24/im flux relation 
for the CHANG-ES galaxies with MIPS images, 
and its obvious correlation further encourages this 
approach. 

To measure 22-/im galaxy isophotal diameters, 
we formed a major axis profile averaged over a 
13" minor axis extent, and measured the diameter 
where the flux had dropped to three times the typ¬ 


ical uncertainty level in the images. This flux level 
is 15.7 /xJy/arcsec^. A measure of the star forma¬ 
tion rate surface density was then calculated by 
assuming disk axisymmetry and dividing the SFR 
by the isophotal disk area. Flux densities, angu¬ 
lar diameters, SFRs and SFR surface densities are 
listed in Table [H 


5.2.1. AGN contamination 


About ten galaxies have bright nuclei at 22 /xm, 
raising the possibility of significant AGN contami¬ 
nation. To give an idea of the potential contamina¬ 
tion, we refer to the study of the AGN contribution 
to mid-IR continuum emission based on observa¬ 


tions of the nuclei of nearby AGN by Tommasin 
et al.j (2010|. They find lower limits of 45-73% for 
the AGN contribution to the 19 /xm continuum, 
depending on the type of AGN. For cases where 
we believe an AGN to be present, we crudely cor¬ 
rect our SFRs by assuming that 100% of the nu¬ 
clear 22/im emission arises from the AGN. A lit¬ 
erature search suggests that all ten candidates are 
starburst nuclei except for two, NGG 3735 and 
NGC 4388, with perhaps the most compelling evi¬ 
dence coming from the IR line ratio study of nuclei 


by Pereira-Santaella et al. (2010). 


We explored this issue further by making a 
rough determination of the nuclear 3.4-4.6 /xm 
colour from the WISE WERGA images, a mea¬ 
sure which distinguishes well between AGN and 
starbursts (Jarrett et al. 2011 Fig. 26, and ref¬ 


erences therein). We estimated the colour of the 
central source in two different ways, both of which 
include simplifying approximations. The resolu¬ 
tions are 5.9" and 6.5" at 3.4 and 4.6 /xm respec¬ 
tively. First, we fit a 2-D Gaussian with a constant 
background to the central peak in each band. We 
tried this with two different box sizes: 9" x 9" 
and 12" xl2". We then turned the flux density 
ratio of the Gaussians into a colour. The second 
method was to Gaussian-convolve the 3.4 /xm im¬ 
ages to the 4.6 /xm resolution, make a colour map, 
and examine the colour in a central box of size 6" 
X 6". The reason these methods are approximate 
is first that the PSFs have significant Airy disks, 
so some of the flux density is not in the central 
peak and not described by a Gaussian, and sec¬ 
ond that a constant background is not necessarily 
the best assumption in the first method since the 
galaxy has structure as it cuts through the fitting 


14 

































box. However, the results are reasonably consis¬ 
tent from the two methods. 


Although Fig. 26 of Jarrett et al. (2011) rep¬ 
resents a range of redshifts, our resulting colours 
clearly fall into either the spiral/starburst or AGN 
regions, and the classifications are generally con¬ 
sistent with the evidence from the literature, with 
two exceptions. 


The first exception is NGC 3735, which has a 
starburst-like nuclear colour of about 0 . 0 - 0 .3 mag 
(the range is from the two different methods and 
shows the greatest discrepancy for any of the ten 


galaxies), whereas MIR line ratios from Pereira- 


Santaella et al. (2010) clearly indicate an AGN. 


In this case, we tend to favour the MIR over the 
NIR evidence because of the possibility of extinc¬ 
tion hiding the NIR emission from the AGN in an 
edge-on geometry. If it is an AGN, we can only es¬ 
timate the SFR by measuring the disk flux density 
outside a radius of 13". This reduces the SFR to 
l.lM 0 /yr, although strictly this is a lower limit. 


The second exception is NGC 660, where a 
colour of 1.2-1.3 clearly indicates an AGN. If in¬ 
deed an AGN, the nuclear flux density can be 
crudely subtracted via a Gaussian ht to the cen¬ 
tral peak in the 22/im image, leaving a lower limit 
SFR of 0.61-0.81 Mq/jt, depending on the fit pa- 


rameters. Yet the literature (e.g. 

Pereira-Santaella 

et al .|2010 

Bernard-Salas et al.||2009 

1 , indicates a 


starburst, as well as a LINER galaxy. Thus, we 
will continue assuming NGC 660 to be a starburst 
rather than AGN galaxy, and retain the originally 
derived values in Table IH 


The nucleus of NGC 4388 does appear to be 
AGN-dominated (see also e.g. Falcke et al.||1998 l, 
with a colour of 1.1 mag. The 22-fim flux density 
is dominated by this source, and the faint disk 
emission is not well resolved from it. Hence, to 
estimate the SFR, we turn to the Spitzer 24-/xm 
image, where the disk is better resolved. However, 
in this case the nuclear source so overwhelms the 
disk that the uncertainty in the flux density of a 
Gaussian fit precludes an estimate of a disk SFR 
as was done for NGC 660. At best, we can crudely 
estimate a disk flux density outside a radius of 7" 
from the nucleus. Hence, the SFR of 0.07 Mq/jy 
is a lower limit. 


5.3. Overview of the CHANG-ES sample 
ordered by SFR surface density 


It has long been suspected that the extent of 
radio halos is related to SFR or SFR per unit area 
(SFR surface density) (e.g. |Dahlem et al. 2006). 
As a preliminary ‘quick-look’ to see whether such 
ordering applies to the CHANG-ES sample, we 
have formed two maps displayed in Figures 
and These maps have been formed from the 
L-band primary-beam-corrected imagej^ They 
have then been ordered by SFR surface density as 
described in the previous section, the latter quan¬ 
tity importantly being independent of distance. 


It is important to keep in mind in this section, 
that in spite of the well-defined criteria used in the 
previous section, some starburst-classified galaxies 
may still have, in addition, hidden AGN contam¬ 
ination, and future examination of the CHANG¬ 
ES sample including higher resolution data sets 
should help to clarify this issue. 

The first map (Fig.|^ is an attempt to bring all 
galaxies to a common distance of 10 Mpc so that 
their 1.6 GHz sizes can be compared. Each image 
was first rotated according to the Ks position an¬ 
gle provided in NED and then minor adjustments 
were made of order < |4°| to produce better hori¬ 
zontal alignment. The galaxies were then scaled in 
size by the ratio, (D/(10 Mpc))^. Appropriate ad¬ 
justments were made to ensure that the brightness 
(i.e. specific intensity) did not change during this 
process. The displayed beams (lower right hand 
corner of each panel) were scaled in the same way 
and, since the resolutions were all similar (though 
not identical) for the L-band data, a glance at 
these beams reveals the distances and linear spa¬ 
tial resolutions of the galaxies; the more distant 
galaxies have large beams after scaling, whereas 
the closer galaxies have small beams. 

Fig. a reveals a wide range of physical sizes for 
the radio continuum emission in the various galax¬ 
ies, both along the disk and into the halo. Some of 
the galaxies show a strong dominant central radio 
core (e.g. NGC 4594), others show a fairly large 
vertical radio extent (e.g. NGC 4666); some ap¬ 
pear rather long and thin (e.g. NGC 4565), others 
show peculiar structure (e.g. NGC 4438 which has 


^®The scaling described in this section was carried out us¬ 
ing the Astronomical Image Processing System (AIPS) of 


NR AO. 
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a radio lobe) and yet others are clearly influenced 
by nearby companions (e.g. NGC 4302). 

The physical resolution, of course, varies for 
each galaxy and a detailed analysis must await 
a future paper. Nevertheless, it is fairly obvious 
that the first map (Fig.|^ is more of a comparison 
of radio continuum sizes than a revelation about 
correlations with SFR per unit area. We also note 
that there is sometimes an uncomfortable range of 
distances in the literature for some of the galaxies. 


The second map (Fig. is again ordered by 
SFR surface density. Flowever, rather than scaling 
by distance, the galaxies are scaled by the 22 /rm 
WISE WERGA major axis sizes (Table |^. The 
galaxy which has the largest angular size on the 
sky is NGC 4244 and the smallest is NGC 5084 
(11.5' and 0.73' respectively). All galaxies were 
therefore scaled up to match the angular size of 
NGC 4244, again ensuring that the brightness was 
not changed and that the beams were also scaled. 
Such a map now essentially corrects for the various 
physical sizes of the galaxies (at 22 /rm). 

In this map, the radio emission still shows some 
size variation because the radio emission embed¬ 
ded within the galaxy might be compact or ex¬ 
tended and, again, the spatial resolution varies. 
Nevertheless, a trend is now apparent, with galax¬ 
ies at the top left (high SFR surface density) 
appearing more ‘boxy’, on average, than galax¬ 
ies at the bottom right (low SFR surface den¬ 
sity). For example, NGC 4244, which has the low¬ 
est SFR/unit area is clearly a very ‘flat’ galaxy 
with no radio halo evident, consistent with the 
H I distribution observed with a deep observation 
by ( Zschaechner et ^|2011 ). On the other hand, 
NGC 4666 which is unlikely to be contaminated 
by an AGN and has good spatial resolution, has 
an obvious halo. 


Explorations of correlations with other param¬ 
eters awaits a future paper. 


6. Median edge-on galaxy and its halo 

The fact that we now have all maps scaled to 
the same 22fim diameter in Fig. allows us to 
form the median edge-on galaxy at L-band. That 
is, we can explore what the typical edge-on spi¬ 
ral galaxy looks like in the radio continuum. We 
have done this for 30 of our 35 galaxies, exclud¬ 
ing NGC 660 (too distorted since it is a merger). 


NGC 4438 (presence of plume/bubble), NGC 4594 
and NGC 5084 (large beam), and UGC 10288 
(emission dominated by a background AGN). 

We formed the median of the 30 galaxies af¬ 
ter having converted the brightness to units of 
Jy/pixel. The result is shown in Fig. and takes 
advantage of the collective sensitivity of all 30 
galaxies. In the image, we also plot a sample 22/im 
contour (red) that corresponds to the scaling of the 
radio data. 


Fig .[^ shows the spectacular radio extent of the 
typical spiral galaxy (note that the average beam 
is much smaller than the displayed radio extent), 
predicted half a century ago by [Ginzburg fc Sy-| 
rovatsky (1961) (see their Fig. 1, p. 18). 


The 

galaxy has the appearance of a slightly flattened 
ellipsoid and reveals that cosmic rays and mag¬ 
netic fields not only permeate the galaxy disk it¬ 
self, but extend far above and below the disk, as 
has been discussed by e.g. [Haverkorn &: Heesenj 
(2012); Krause (2009). We note that the red el¬ 
lipse is not affected by the larger beams of the ra¬ 
dio images and should not be used, at this point, 
to make conclusions about the radial distribution 
of the emission compared to the 22/im emission. 
On the other hand, the conclusion of a broadscale 
halo is robust; we have checked the contribution 
of the disk emission to the apparent vertical radio 
extent by examining the lower inclination galaxies 
in the sample, finding that the disk made a negligi¬ 
ble contribution. If galaxies that have lower incli¬ 
nation are removed, the halo remains, suggesting 
that the vertical distribution is not a result of a 
projected disk. In addition, the halo extends be¬ 
yond what would be expected from beam smooth¬ 
ing. 


The galaxy scale heights (see Dumke et al. 1995 


for method) and magnetic field extents will be dis¬ 
cussed further in future papers. 


7. Conclusions 

1. In this CHANG-ES IV paper, we present all 
the VLA D-configuration observations and 
results of the 35 edge-on CHANG-ES galax¬ 
ies in two frequency bands, C and L. These 
data products (including intensity maps, 
spectral index maps, polarization maps) are- 
all be part of our Data Release 1, located 
at http://www.queensu.ca/changes. Apart 
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from presenting each galaxy in a range of 
maps for two different weightings of each 
band, we also investigate deeper the in-band 
spectral index maps, primary beam correc¬ 
tions, flux densities and star formation rates. 

2. Spectral index maps presented have a) been 
primary beam corrected, b) cut off below 
5(7, c) cut off wherever the formal error > 
1.0, and been convolved with a Gaussian 
to smooth over the pixel-to-pixel variations 
within any clean beam. No changes in spa¬ 
tial resolution results from this process. 

3. Spectral index uncertainties have been thor¬ 
oughly investigated. Apart from extreme 
values near the edges ’’edge-effect” which 
should be ignored, there are effects of the 
primary beam corrections. Our tests on two 
pointing data show that results outside the 
half-power level of the primary beam model 
used in CASA need be treated with cau¬ 
tion. The true error in a is typically under¬ 
estimated by 20% in the formal error maps. 
In regions of high signal-to-noise, they may 
even underestimate the error by a factor ~5. 
The latter affects C-band more than L-band, 
since the galaxies do not extend beyond the 
half-power point of the primary beam for L- 
band. 

4. Galaxies with flatter a values at the centres 
often have indications of central activity. In 
some cases, such as NGC 2820, NGC 3556 
and NGC 5775 and others, a is flatter at 
C-band than in L-band, suggesting a higher 
contribution from thermal emission at the 
higher frequency. 

5. Star formation rates were derived via WISE- 
WERGA images, and potential AGN con¬ 
tamination of the SFR results investigated. 
Many candidates may have starburst nuclei, 
with the exceptions of NGC 660, NGC 3735 
and NGC 4388 who may indeed harbour an 
AGN. Lower limit SFR are calculated for 
these and used for the two latter galaxies. 

6. We scaled, rotated and ordered the L-band 
images of the CHANG-ES galaxies by SFR 
surface density in order to be able to com¬ 
pare the galaxies on equal footing and get a 


snap shot of potential correlations between 
SFR and halo size. We produced two maps, 
one scaled by distance (to 10 Mpc), the other 
by the 22^m WISE-WERGA sizes, and the 
second one does reveal a trend of galaxies 
harbouring larger halos having higher SFR 
densities. 

7. The L-band images from the previous point 
were medianed to bring forth the median 
edge-on galaxy, which harbours a compelling 
halo. 
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A. A smorgasbord of results for each CHANG-ES galaxy 


The figures in this appendix display an assortment of the data products from the CHANG-ES project, 
described below. With the exception of the field panel (j) in the red frame, all panels are displayed in the 
exact same angular field of view, for better comparison. The ’apparent B-vectors’ in panels (d), (e) and (f) 
below, are in reality E-vectors rotated by 90 degrees (see Sec. 3.3.3). 


First row 

a) C-band Stokes I intensity map with a Briggs robust 0 weighting. The I contours are overlaid for clarity, 
at 3, 6, 12, 24, 48, 96, 192 times the rms noise of the map (rms values are listed Table [^. 

b) C-band Stokes I intensity map with a uv-tapering applied onto the robust 0 weighting, yielding a beam 
approximately twice the size of that of the robust 0 weighting in C-band, overlaid with its contours at 3, 6, 
12, 24, 48, 96, 192 times the rms noise. This map does not exist for galaxies NGC 4438 and NGC 4845. 

c) L-band Stokes I intensity map, with contours at 3, 6, 12, 24, 48, 96, 192 times the L rms noise (see Table[^ 
for rms values). 


Second row: 

d) C-band Stokes I contours with uv-tapered weighting at 3 times the rms noise, with apparent B-vectors at 
the same weighting , overlaid on an optical image. The optical images were created using a combination of 
Sloan Digital Sky Survey (SDSS) g,r,i bands or Digitized Sky Survey 2 (DSS2) blue red and infrared bands, 
for the galaxies not available in SDSS. 

e) C-band polarization map with apparent B-vectors and Stokes I contours at 3 ct. 

f) L-band polarization map with apparent B-vectors and Stokes I contours at 3cr. 


Third row: 

g) Stokes I intensity map with 3 weightings merged (images a, b and c above). Image a, i.e. C-band with 
robust 0 weighting, is shown in dark blue, image (b), i.e. C-band with uv-tapering, is shown in green, and 
image c, L-band robust 0 weighting, is shown in white. Except for six galaxies for which no satisfactory 
uv-tapered L-band image has been obtained (see Table the background has been masked out using the 
L-band uv-tapered map (this map is not shown separately). 

h) C-band spectral index map. Spectral indices have been cut off where formal errors are larger than 1.0, 


and then smoothed/averaged to the size of the beam (Section 3.4). The colour bar of panel i) also applies 
to this panel. 

i) L-band spectral index map. Spectral indices have been cut off where formal errors are larger than 1.0, 
and then smoothed/averaged to the size of the beam (Section |3.4[). 


Fourth row: 

j) The L-band robust 0 surrounding environment. The red frame indicates that this panel is displayed in a 
different field of view than the other panels. 

k) Spectral index uncertainty map of (h) - formal error data. See Section 3.4.3 for more information on the 


uncertainties. Note that although the true error is higher than what is displayed in this map, the map gives 
an indication of which regions are to be treated with caution. The colour bar of panel (1) also applies to this 
panel. 


1) Spectral index uncertainty map of (h) - formal error data. See caption of panel (k) and Section 3.4.3 for 
more information on the uncertainties. 
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Fig. 1.— a) Absolute difference between the PB-corrected spectral index measurements in the two C-band 
pointings of NGC 891. There is an increase in difference with distance from the midpoint between the 
pointing centres, even in regions with high S/N. b) Error of the PB-weighted average of the PB-corrected 
spectral index measurements in the two pointings. The red contour is placed at an error of 0.1. c) Spectral 
index error map determined by the ms-mfs clean algorithm (same as panel (k) of the appendix figure of 
NGC 891, but without the cutoff at Aa > 1.0). Contours represent an error of 0.1 in this map (white) and 
in the map shown in panel b) (red), d) Ratio of the two different error maps (panel b) divided by panel c)). 
The ratio is mostly close to unity, except for those parts of the disk lying outside the 0.7 PB level of either 
pointing. In each panel, the two PB circular contours for each pointing are placed at the 0.5 and 0.7 level. 
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Fig. 2.— Construction of combined weightings and bands in panel (g) of the figures in Appendix. The layering 
procedure in the Gnu Image Manipulation Program (GIMP) is used to combine the available weightings of 
the C and L band data. (”robO” stands for the robust 0 weighting and ’’uvtap” is the uv-tapered weighting.) 
The algorithmic mode applied to the top layer is listed under the upper image in each column. The resulting 
combination is displayed in the bottom row. The left hand column represents the first 2 layers that are 
combined, i.e. inverted-intensity and colourized versions of the two weightings of the C band data. The 
second step, represented in the middle column, combines the inverted-intensity higher resolution L band 
data with the result of Step 1. The right hand column shows how the inverted-intensity, lower resolution 
L band data are applied to the result of Step 2 in order to mask out confusing structures in the off-target 
’’background”. The bottom image in the right hand column is presented as panel (g). 
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Fig. 4.— The CHANG-ES galaxies L-band data, rotated and scaled by size to a distance of 10 MPc, ordered 
by SFR density. Highest SFR density at top left, decreasing by row, left to right. Contours shown are at 
levels of 1, 5, 25 and 125 mJy/beam. 5' corresponds to 14.5 kpc at this scale. 
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Fig. 5.— The CHANG-ES galaxies L-band data, rotated, and scaled to the 22 micron WISE-WERGA size 
of NGC 4244, ordered by SFR surface density. Highest SFR surface density at top left, decreasing by row, 
left to right. Contours shown are at levels of 1, 5, 25 and 125 mJy/beam. 
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Fig. 6.— The median edge-on spiral galaxy in L-band, made from stacking 30 of the galaxies in Fig. The 
red ellipse is a sample 22/im contour that corresponds to the scaling of the radio data and thus represents 
the disk radial extent. The beam shown is the average beam of the 30 galaxies. 
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Galaxy®" RA dec Distance^ (Mpc) Band Date (yymmdd) SB ID prim. cal. ^ zero. pol. cal.^ sec. cal. 

N 660 01h43m02.s40 +13d38m42.s2 12.3 C 111209 5804098 3C48 3C84 J0204+1514 

C 111218 5805266 3C48 3C84 J0204+1514 
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L 111230 4812476 3C286 OQ208 J1254+1141 

N 4217 12hl5m50.s90 +47d05m30.s4 20.6 C 111227 5062561 3C286 OQ208 J1219+4829 
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Table 2 

Two POINTINGS OF LARGE GALAXIES 


Galaxy 

(1) 

RA 1 
(2) 

Dec 1 

(3) 

RA 2 

(4) 

Dec 2 

(5) 

N891 

02h 22m 37.21s 

42d 22m 41.2s 

02h 22m 29.61s 

42d 19m 12.6s 

N3628 

llh 20m 24.50s 

13d 34m 55.7s 

llh 20m 9.52s 

13d 35m 50.1s 

N4244 

12h 17m 36.71s 

37d 49m 40.9s 

12h 17m 22.61s 

37d 47m 10.3s 

N4565 

12h 36m 26.58s 

25d 57m 54.7s 

12h 36m 14.98s 

26d 00m 36.6s 

N4631 

12h 42m 16.88s 

32d 32m 37.2s 

12h 41m 59.14s 

32d 32m 21.6s 

N4594 

12h 40m 07.09s 

- lid 37m 23.0s 

12h 39m 51.77s 

- lid 37m 23.0s 

N5084 

13h 20m 24.88s 

- 21d 49m 19.8s 

13h 20m 8.96s 

- 21d 49m 58.8s 

N5907 

15h 15m 59.49s 

56d 18m 1.6s 

15h 15m 48.05 

56d 21m 25.5s 


Note.— Column 1: Galaxy name; Column 2: Right Ascension of pointing 1; Column 
3: Declination of pointing 1; Column 4: Right Ascension of pointing 2; Column 5: 
Declination of pointing 2. 


Table 3 


Peeled sources in the 


L-band fields of two galaxies 


Galaxy 

Source name 

Angular distance 

Intensity before/after 

RA 

Dec 



(') 

(mjy/beam) 



(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

N660 

_ 

15 

95/26 

Olh 42m 18.8s 

+ 13d 27m 48.6s 

N4388 

M87 

75 

132/22 

12h 30m 48s 

+ 12d 23m 30.2s 


Note.— Column 1: Galaxy name; Column 2: Name of interfering source 1; Column 3: distance in arcminutes 
between galaxy and interfering source; Column 4: Intensity of interfering source before/after peeling ; Column 
5: Right Ascension; Column 6: Declination. 
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Table 4 

Imaging results C-band 


Galaxy 

(1) 

weighting 

( 2 ) 

Bmaj 

n 

(3) 

^min 

[”] 

( 4 ) 

BPA 

[degree] 

(5) 

Vo 

[GHz] 

(6) 

rms I 

[/ijy/beam] 

(7) 

Map peak 
[mjy/beam] 
(8) 

DR 

(9) 

rms QU 
[^Jy/beam] 
(10) 

pol. peak 
[/j.Jy/beam] 
(11) 

N660 

rob 0 

9.55 

9.14 

-18.72 

5.99904 

8.6 

603 

70116 

6.0 

369.1 


uvtap 

16.13 

16 

-58.32 


18 

614 

34111 

7.0 

554.6 

N891 

rob 0 

9 

8.81 

-79.46 

5.99932 

6.5 

7.21 

1109 

6.6 

115.0 


uvtap 

15.31 

14.63 

88.05 


6.9 

6.9 

1000 

7.2 

225.0 

N2613 

rob 0 

15.42 

8.3 

-11.05 

5.99900 

10.1 

0.643 

64 

10.4 

57.1 


uvtap 

19.94 

15.82 

-11.44 


11.8 

1.09 

92 

11.5 

63.6 

N2683 

rob 0 

9.36 

8.76 

11.71 

5.99900 

9.1 

3.27 

359 

8.8 

50.2 


uvtap 

15.88 

14.81 

46.67 


9.4 

3.31 

352 

8.9 

64.6 

N2820 

rob 0 

9.84 

8.91 

-3.81 

5.99900 

5.3 

1.699 

321 

5.3 

48.7 


uvtap 

16.18 

15.36 

-3.884 


5.2 

3.173 

610 

5.3 

97.3 

N2992 

rob 0 

14.33 

8.84 

-13.3 

5.99900 

8.4 

65.6 

7810 

8.0 

151.7 


uvtap 

18.51 

16.76 

-14.26 


9.5 

71.7 

7547 

8.2 

202.8 

N3003 

rob 0 

9.37 

8.93 

-67.02 

5.99900 

9.4 

1.35 

144 

8.6 

41.1 


uvtap 

13.78 

13.61 

16.65 


6.0 

8.11 

1352 

9.2 

40.0 

N3044 

rob 0 

10.68 

9.59 

-5.08 

5.99900 

5.1 

5.92 

1161 

7.0 

102.9 


uvtap 

13.78 

13.61 

16.65 


6.0 

8.11 

1352 

7.0 

138.4 

N3079 

rob 0 

9.32 

8.72 

-13.89 

5.99900 

6.5 

202.07 

31088 

5.9 

1987.1 


uvtap 

16.16 

15.44 

-9.5 


8.0 

221.93 

27741 

5.8 

2759.5 

N3432 

rob 0 

11.94 

9.97 

3.43 

5.99900 

7.2 

2.33 

324 

7.1 

39.2 


uvtap 

15 

13.14 

-179.3 


6.8 

2.62 

385 

6.5 

49.3 

N3448 

rob 0 

9.1 

8.6 

3.4 

5.99900 

5.6 

5.36 

957 

5.4 

66.0 


uvtap 

16 

15.7 

-167.6 


6.6 

8.46 

1282 

5.9 

105.5 

N3556 

rob 0 

9.2 

8.7 

7.6 

5.99900 

5.3 

3.45 

651 

5.3 

66.8 


uvtap 

16 

15. 

14.1 


6.0 

6.9 

1150 

5.7 

91.9 

N3628 

rob 0 

9.62 

9.47 

-31.09 

5.99842 

9.5 

66.96 

7048 

10.2 

174.2 


uvtap 

16.83 

15.15 

75.57 


9.5 

79.09 

8325 

9.6 

128.4 

N3735 

rob 0 

10.1 

8.8 

-2.3 

5.99900 

5.8 

2.9 

500 

5.6 

109.5 


uvtap 

16.4 

15.6 

-5 


6.2 

5.3 

855 

5.6 

220.4 

N3877 

rob 0 

8.85 

8.66 

39.45 

5.99900 

6.5 

1.43 

220 

6.4 

25.6 


uvtap 

16.07 

15.27 

22.61 


6.7 

1.94 

290 

6.6 

33.7 

N4013 

rob 0 

9.07 

8.88 

-31.87 

5.99900 

6.4 

3.88 

606 

6.4 

41.5 


uvtap 

16.06 

15.25 

-7.25 


6.6 

4.58 

694 

6.5 

62.5 

N4096 

rob 0 

9 

8.85 

-59.54 

5.99900 

6.3 

0.529 

84 

6.3 

26.9 


uvtap 

15.85 

15.32 

-16.41 


6.8 

1.1 

162 

6.6 

28.7 

N4157 

rob 0 

9.1 

8.74 

44.64 

5.99900 

6.3 

1.47 

3411 

6.0 

58.3 


uvtap 

15.9 

14.99 

21.61 


6.2 

3.6 

581 

6.1 

123.3 

N4192 

rob 0 

9.17 

8.99 

-11.55 

5.99900 

6.0 

3.807 

635 

6.0 

150.2 


uvtap 

15.66 

15 

83.71 


6.0 

5.128 

855 

5.9 

222.2 

N4217 

rob 0 

8.91 

8.73 

29.41 

5.99900 

6.3 

2.63 

417 

6.3 

69.4 


uvtap 

16.06 

15.23 

17.5 


6.3 

4.52 

717 

6.2 

137.5 

N4244 

rob 0 

9.25 

8.91 

-4.37 

5.99838 

5.6 

0.693 

235 

5.9 

14.3 


uvtap 

15.83 

15.07 

-1.57 


5.8 

0.84/1.39 

240 

5.7 

18.6 

N4302 

rob 0 

9.96 

8.96 

-21.41 

5.99900 

15.5 

1.29 

146 

15.2 

73.1 


uvtap 

18.7 

18.4 

32.98 


19.5 

1.93 

99 

17.0 

113.3 

N4388 

rob 0 

9.9 

8.99 

-24.75 

5.99900 

13.9 

18.1 

1302 

13.5 

436.1 


uvtap 

16.2 

15.97 

-23.7 


18.0 

25.17 

1398 

16.4 

462.7 

N4438 

rob 0 

9.83 

9.17 

-19.86 

5.99900 

16.0 

22.91 

1432 

15.0 

221.8 

N4565 

uvtap n/a 
rob 0 

9.02 

8.82 

85.48 

5.99834 

7.4 

2.64 

357 

7.4 

49.7 


uvtap 

15.96 

14.83 

42.18 


8.0 

3.11 

389 

7.8 

105.0 

N4594 

rob 0 

13.32 

8.91 

1.1 

5.99841 

12.9 

102 

7907 

11.4 

228.2 


uvtap 

18.09 

16.07 

0.8 


14.4 

103.3 

7159 

13.7 

245.0 

N4631 

rob 0 

8.88 

8.57 

-6.88 

5.99835 

7.7 

6.8 

883 

6.9 

122.4 


uvtap 

15.71 

14.73 

86.79 


9.2 

13.9 

1511 

9.6 

368.6 

N4666 

rob 0 

11.06 

9.5 

-179.82 

5.99900 

12.5 

7 

560 

12.0 

210.8 
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Table 4 —Continued 


Galaxy 

(1) 

weighting 

(2) 

Bmaj 

n 

(3) 

Bjnin 

[”] 

(4) 

BPA 

[degree] 

(5) 

1^0 

[GHz] 

(6) 

rms I 

[/ijy/beam] 

(7) 

Map peak 
[mjy/beam] 
(8) 

DR 

(9) 

rms QU 
[^Jy/beam] 
(10) 

pol. peak 
[/j.Jy/beam] 
(11) 


uvtap 

14.1 

13.54 

-176.99 


12.5 

10.61 

849 

11.5 

332.0 

N4S45 

rob 0 

10.98 

9.06 

-1.4 

5.99900 

15.0 

424.7 

28313 

14.5 

337.3 


uvtap n/a 










N5084 

rob 0 

15.64 

8.35 

-5.76 

5.99841 

7.9 

29.1 

3684 

8.0 

76.4 


uvtap 

20.33 

16.23 

-11.18 


8.6 

29.4 

3419 

9.4 

83.3 

N5297 

rob 0 

9.03 

8.78 

35.28 

5.99900 

6.3 

0.267 

42 

5.6 

24.9 


uvtap 

15.66 

14.87 

21.77 


6.1 

0.485 

SO 

6.1 

29.8 

N5775 

rob 0 

10.08 

9.34 

-22.79 

5.99900 

5.0 

4.45 

890 

5.0 

104.0 


uvtap 

16.1 

14.96 

75.53 


5.8 

7.77 

1340 

5.3 

214.2 

N5792 

rob 0 

10.34 

9.18 

-12.72 

5.99900 

5.3 

6.183 

1167 

5.3 

73.9 


uvtap 

16.09 

15.54 

-87.29 


5.4 

9.29 

1720 

5.4 

76.9 

N5907 

rob 0 

9.78 

8.08 

16.32 

5.99847 

9.3 

0.875/5.76 

623 

9.1 

140.7 


uvtap 

16.65 

14.55 

65.68 


9.8 

1.17/5.22 

533 

9.6 

115.5 

U10288^ 

^ rob 0 

10.96 

9.46 

-28.47 

5.99900 

7.0 

8.94 

1277 

5.4 

146.6 


uvtap 

13.99 

13.42 

-59.75 


8.5 

9.39 

1105 

5.4 

169.8 

Note.— 

-Col. 1: Galaxy 

name; 

Col. 2: Weighting, where ”rob 0” 

refers to Briggs weighting with robust value set to 0, 

’’uvtap” refers 


to a uv-tapered version of the former; Columns 3-5: synthesized beam parameters; Col. 6: C-band central frequency in GHz (varying 
due to differences in flagging); Col. 7: Stokes I rms noise; Col. 8: Peak intensity of the galaxy. In the case of two values given, the first 
value is the peak intensity of the galaxy and the second higher value is the peak of the map when this occurs outside of the galaxy; Col. 
9: Dynamic range in image (map peak intensity over noise); Col. 10: Stokes Q and U average rms noise; Col. 11: Peak intensity of the 
polarization map (measured from non-PB-corrected maps). 

^Map peak and polarization peak values given for UGC 10288 are of the background source. 


66 



Table 5 

Imaging results L-band 


Galaxy 

(1) 

weighting 

(2) 

Bmaj 

n 

(3) 

[”] 

(4) 

BPA 

[degree] 

(5) 

[GHz] 

(6) 

rms I 

[/ijy/beam] 

(7) 

Map peak 
[mjy/beam] 
(8) 

DR 

(9) 

rms QU 
[^Jy/beam] 
(10) 

pol. peak 
[/j.Jy/beam] 
(11) 

N660 

rob 0 

39.78 

36.27 

80.22 

1.57502 

60 

422.5 

7042 

28 

191 


uvtap 

49.85 

47.29 

-86.89 


75 

441.8 

5891 

28 

177 

N891 

rob 0 

36.42 

32.48 

-74.29 

1.57498 

60 

73.4/78.4 

1307 

40^ 

455 


uvtap 

44.94 

43.11 

67.61 


95 

96.7 

1018 

44 

647 

N2613 

rob 0 

73.43 

41.97 

-29.89 

1.57497 

42 

11.96/17.2 

410 

35 

176 


uvtap 

76.55 

51.18 

-32.11 


45 

13.4/16.9 

376 

40 

191 

N2683 

rob 0 

33.9 

31.31 

-43.8 

1.57481 

30 

8.6 

287 

27 

223 


uvtap 

44.45 

41.69 

-44.23 


35 

11.7 

334 

27 

218 

N2820 

rob 0 

34.8 

30.56 

-32.93 

1.57479 

40 

18.54/39.9 

998 

32 

132 


uvtap 

43.97 

39.88 

-36.87 


40 

23.6 

590 

32 

152 

N2992 

rob 0 

52.08 

33.1 

-13.55 

1.57481 

30 

193.2 

6440 

26 

139 


uvtap 

59.34 

45.81 

-20.87 


32 

195.4 

6106 

27 

167 

N3003 

rob 0 

33.86 

30.7 

-46.95 

1.57478 

30 

8.93 

298 

27 

102 


uvtap 

44.64 

40.58 

-47.84 


29 

11.5 

397 

27 

92 

N3044 

rob 0 

41.79 

32.25 

-48.11 

1.57479 

28 

43.96 

1570 

25 

246 


uvtap 

50.92 

42.39 

-61.13 


27 

20.76 

769 

27 

270 

N3079 

rob 0 

34.88 

31.49 

-48.61 

1.57477 

45 

427.3 

9496 

30 

488 


uvtap 

45.64 

41.61 

-48.18 


45 

481.2 

10693 

30 

646 

N3432 

rob 0 

32.62 

31.82 

-64.03 

1.57474 

36 

12.16/27.83 

773 

32 

181 


uvtap 

42.34 

41.43 

-51.85 


40 

16.8/27 

675 

32 

195 

N3448 

rob 0 

34.55 

31.4 

-23.24 

1.57475 

35 

25.46 

727 

26 

114 


uvtap n/a 










N3556 

rob 0 

34.45 

31.31 

-22.75 

1.57475 

38 

33.87 

891 

25 

372 


uvtap 

45.91 

41.74 

-21.92 


42 

48.2 

1148 

26 

503 

N3628 

rob 0 

32.27 

32.65 

-37.4 

1.57473 

49 

249 

5082 

30 

30 


uvtap 

47.19 

42.47 

-41.85 


55 

267 

4855 

28 

40 

N3735 

rob 0 

36.8 

32.1 

-15.7 

1.57477 

32 

33.5 

1047 

25 

213 


uvtap 

46.6 

41.1 

-21.7 


31 

42 

1355 

25 

250 

N3877 

rob 0 

34.12 

32.61 

-31.6 

1.57473 

30 

7.93 

264 

27 

180 


uvtap 

45.41 

41.36 

-27.75 


37 

10.2 

276 

25 

188 

N4013 

rob 0 

34.15 

31.54 

-65.96 

1.57472 

32 

13.8 

431 

27 

103 


uvtap 

44.98 

41.45 

-64.79 


38 

15.5/24.6 

647 

26 

116 

N4096 

rob 0 

33.85 

31.73 

-35.65 

1.57473 

30 

7.099 

237 

27 

133 


uvtap 

44.75 

41.51 

-34.28 


32 

10.33 

323 

26 

180 

N4157 

rob 0 

34.12 

31.65 

-23.32 

1.57488 

29 

31.42/53.36 

1840 

28 

254 


uvtap 

45.19 

41.48 

-21.07 


35 

43.8/54.5 

1557 

27 

291 

N4192 

rob 0 

35.52 

32.66 

-52.17 

1.57471 

40 

17.5 

438 

28 

205 


uvtap n/a 










N4217 

rob 0 

34.01 

31.46 

-29.04 

1.57472 

28 

24 

857 

27 

264 


uvtap 

45.31 

41.07 

-26.28 


30.4 

33.1 

1089 

28 

368 

N4244 

rob 0 

33.98 

32.46 

-47.1 

1.57471 

30 

2.18/9.87 

329 

27 

114 


uvtap 

45.43 

41.84 

-43.9 


30 

10.1 

337 

27 

93 

N4302 

rob 0 

35.82 

33.55 

-70.51 

1.57470 

46 

9.73/15.32 

333 

35 

319 


uvtap 

47 

42.8 

-81.78 


60 

15.11 

252 

35 

409 

N4388 

rob 0 

36.32 

33.24 

-48.92 

1.59888 

150 

79.92/296.2 

1975 

65 

436 


uvtap n/a 










N4438 

rob 0 

34.7 

30.21 

-53.59 

1.77468 

130 

79.2 

609 

50 

237 


uvtap n/a 










N4565 

rob 0 

34.5 

32.32 

-89.06 

1.57470 

30 

9.36/48 

1600 

27 

270 


uvtap 

43.27 

41.92 

-83.71 


32 

12.58/48.9 

1528 

24 

451 

N4594 

rob 0 

47.89 

32.62 

-4.6 

1.57471 

31 

81.4 

2626 

25 

209 


uvtap 

54.41 

44.46 

-8.23 


32 

81.9 

2559 

25 

254 

N4631 

rob 0 

35 

32.4 

-89.34 

1.57488 

31 

90.1 

2906 

28 

353 


uvtap 

32.81 

41.72 

-82.91 


31 

123 

3955 

33 

449 

N4666 

rob 0 

37.4 

36.08 

-30.42 

1.57470 

23 

102.4 

4452 

25 

700 


uvtap 

47.95 

44.7 

-72.66 


27 

132 

4981 

25 

865 
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Table 5 —Continued 


Galaxy 

(1) 

weighting 

(2) 

Bmaj 

n 

(3) 

Bm.i„ 

[”] 

(4) 

BPA 

[degree] 

(5) 

I'D 

[GHz] 

(6) 

rms I 

[/ijy/beam] 

(7) 

Map peak 
[mjy/beam] 
(8) 

DR 

(9) 

rms QU 
[/j.Jy/beam] 
(10) 

pol. peak 
[//Jy/beam] 
(11) 

N4845 

rob 0 

38.58 

34.27 

-5.22 

1.57470 

40 

224.8 

5620 

27 

556 

N5084 

uvtap 
rob 0 

57.12 

31.48 

-5.85 

1.57472 

33 

35.6 

1079 

29 

395 


uvtap 

65.86 

43.94 

-8.33 


35 

36.5 

1043 

30 

426 

N5297 

rob 0 

34.33 

32.17 

5.23 

1.57471 

34 

4.25 

125 

27 

138 


uvtap 

45.16 

41.38 

-5.49 


32 

6.11 

191 

26 

143 

N5775 

rob 0 

40.35 

35.4 

-42.84 

1.57468 

35 

57.5 

1643 

26 

542 


uvtap 

51.16 

45.04 

-63.19 


33 

75.5 

2323 

27 

645 

N5792 

rob 0 

40.13 

35.23 

-29.48 

1.57469 

32 

31.92 

998 

25 

1366 


uvtap 

49.54 

47.71 

-39.27 


35 

34.4 

983 

25 

1461 

N5907 

rob 0 

32.99 

30.67 

44.87 

1.57474 

30 

9.45/52.79 

1759 

28 

556 


uvtap 

42.49 

39.93 

42.4 


30 

13.4/31.8 

1060 

28 

603 

U10288 

rob 0 

40.23 

34.34 

-31.32 

1.57488 

39 

98.7 

2531 

35 

3027 


uvtap n/a 










Note.- 

-Col. 1: Galaxy 

name; 

Col. 2: Weighting, where ”rob 0” 

denotes Briggs weighting with robust value set to 0, 

’’uvtap” refers 


to a uv-tapered version of the former; Columns 3-5: synthesized beam parameters; Col. 6: L-band central frequency in GHz (varying 
due to differences in flagging); Col. 7: Stokes I rms noise; Col. 8: Peak intensity of the galaxy. In the case of two values given, the first 
value is the peak intensity of the galaxy and the second higher value is the peak of the map; Col. 9: Dynamic range in image; Col. 10: 
Stokes Q and U average rms noise; Col. 11: Peak intensity of the polarization map (measured from non-PB-corrected maps). 

previous observation (an extra 10 minutes on source) was included in the L-band Stokes Q and U imaging to increase sensitivity. 
It was however excluded from Stokes I imaging due to artifact contamination. 

*^Map peak and polarization peak values given for UGC 10288 are of the background source. 
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Table 6 

Star formation rates 


Galaxy 

(1) 

Flux density 

[Jy] 

(2) 

Flux 

[xW-^^W/rrr‘\ 

(3) 

Diameter 

['] 

(4) 

Diameter 

[kpc] 

(5) 

SFR 
[Mq/ yr] 
(6) 

SFR surface density 
[x 10~^Mq {yrjkpc^] 

(T) 

Uncertainty 

(8) 

N660 

5.70 

7.77 

3.02 

10.8 

2.74 

29.96 

0.02 

N891 

5.90 

8.04 

9.50 

25.1 

1.55 

3.13 

0.02 

N2613 

0.99 

1.35 

4.95 

33.7 

1.73 

1.94 

0.02 

N26S3 

0.74 

1.01 

5.15 

9.39 

0.09 

1.34 

0.02 

N2820 

0.28 

3.80 

1.78 

13.8 

0.62 

4.20 

0.02 

N2992 

0.88 

1.19 

1.27 

12.5 

3.22 

26.13 

0.03 

N3003 

0.33 

4.46 

3.75 

27.7 

0.67 

1.11 

0.02 

N3044 

0.73 

9.91 

3.07 

18.1 

0.95 

3.70 

0.02 

N3079 

2.56 

3.50 

4.33 

26.0 

3.46 

6.54 

0.02 

N3432 

0.54 

7.42 

3.63 

9.96 

0.15 

1.97 

0.02 

N3448 

0.48 

6.59 

1.75 

12.5 

0.92 

7.55 

0.02 

N3556 

3.43 

4.68 

6.08 

24.9 

2.17 

4.44 

0.02 

N3628 

4.41 

6.02 

9.90 

24.5 

1.01 

2.15 

0.02 

N3735 

0.84 

1.14 

2.82 

34.4 

1.10^ 

1.2^ 

0.02 

N3877 

0.92 

1.25 

3.58 

18.5 

0.92 

3.43 

0.02 

N4013 

0.59 

7.98 

3.45 

16.1 

0.48 

2.35 

0.02 

N4096 

0.79 

1.08 

3.92 

11.8 

0.27 

2.46 

0.02 

N4157 

1.62 

2.21 

3.67 

16.6 

1.25 

5.77 

0.02 

N4192 

0.96 

1.31 

6.20 

24.4 

0.56 

1.19 

0.02 

N4217 

1.14 

1.55 

3.90 

23.4 

1.53 

3.57 

0.02 

N4244 

0.38 

5.22 

11.53 

14.8 

0.02 

0.14 

0.03 

N4302 

0.44 

6.00 

3.82 

21.6 

0.53 

1.45 

0.02 

N4388 

2.18 

2.97 

2.30 

11.1 

0.07*= 

0.72*= 

0.02 

N4438 

0.19 

2.64 

1.32 

3.98 

0.07 

5.35 

0.03 

N4565 

1.65 

2.25 

10.43 

36.1 

0.74 

0.73 

0.02 

N4594 

0.62 

8.44 

5.97 

22.0 

0.32 

0.83 

0.02 

N4631 

7.62 

1.04 

10.85 

23.4 

1.33 

3.10 

0.02 

N4666 

3.03 

4.13 

4.03 

32.3 

7.29 

8.92 

0.02 

N4845 

0.52 

7.14 

2.10 

10.4 

0.48 

5.68 

0.02 

N50S4 

0.06 

7.73 

0.73 

4.99 

0.10 

5.04 

0.05 

N5297 

0.24 

3.33 

2.20 

25.9 

1.27 

2.41 

0.02 

N5775 

1.99 

2.71 

3.80 

32.0 

5.28 

6.58 

0.02 

N5792 

0.82 

1.12 

2.57 

23.7 

2.63 

5.98 

0.02 

N5907 

1.73 

2.36 

7.35 

35.9 

1.56 

1.54 

0.02 

U10288 

0.11 

1.51 

2.15 

21.3 

0.41 

1.15 

0.03 


Note.— Col. 1: Name; Col. 2 and 3: 22 fim flux, given in two different units (The values in Col. 3 were used for Fig. ; 
Col. 4 and 5: Angular diameters; Col.6: Star formation rates in solar masses/year; Col 7: Star formation rate density Xl0~^ 
in solar masses/year and kpc^\ Col. 8: Fractional error on surface density. 

^The SFR and SFR surface density values of NGC 3735 were adju sted fr om 4.71 and 0.005 respectively, to the values given 
here (lower limits), in order to account for the central AGN (see Sec. |5.2~T] l 

^The SFR and SFR surface density values of NGC 4388 were adju sted fr om 1.91 and 0.020 respectively, to the values given 
here (lower limits), in order to account for the central AGN (see Sec. |5.2~T|l 
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Table 7 
Flux densities 


Galaxy 

(1) 

Flux density C 

(2) 

Flux density G uncertainty 

(3) 

Flux density L 

(4) 

Flux density L uncertainty 

(5) 

N660 

657.9 

13.2 

525.4 

10.5 

N891 

208.7 

6.9^ 

743.9 

14.9 

N2613 

15.3 

0.3 

59.6 

2 4a 

N2683 

20.3 

0.8^ 

66.6 

6.5^ 

N2820 

19.1 

0.4 

61.8 

1.2 

N2992 

80.4 

1.6 

204.9 

4.1 

N3003 

10.8 

0.5^ 

34.9 

0.7 

N3044 

37.5 

0.9^ 

104.2 

2.1 

N3079 

365.4 

7.3 

811.0 

16.2 

N3432 

26.3 

0.5 

83.3 

1.9^ 

N3448 

20.5 

0.4 

46.0 

0.9 

N3556 

79.2 

4.7^ 

291.5 

5.8 

N3628 

184.6 

3.7 

527.5 

10.5 

N3735 

24.9 

0.5 

81.3 

1.6 

N3877 

12.9 

0.3^ 

42.7 

0.9 

N4013 

12.6 

0.3 

37.9 

0.8^ 

N4096 

16.3 

0.3 

57.1 

1.1 

N4157 

55.1 

1.1 

184.5 

3.7 

N4192 

24.4 

0.5 

80.6 

1.6 

N4217 

35.4 

0.7 

111.5 

2.2 

N4244 

9.0 

0.6^ 

18.1 

0.6^ 

N4302 

12.0 

0.3^ 

45.1 

0.9 

N4388 

62.2 

1.2 

130.9 

2.6 

N4438 

54.5 

1.1 

132.2 

2.6 

N4565 

42.3 

1.1^ 

152.2 

3.0 

N4594 

128.3 

2.6 

93.7 

1.9 

N4631 

284.4 

7.4^ 

1083.0 

37.0^ 

N4666 

125.3 

2.5 

404.5 

8.1 

N4845 

432 

8.6 

230.0 

4.6 

N5084 

36.1 

1.0^ 

40.7 

0.8 

N5297 

6.7 

0.4^ 

24.4 

0.5 

N5775 

74.4 

1.5 

255.0 

5.1 

N5792 

20.2 

0.4 

57.7 

1.2 

N5907 

51.5 

1.0 

181.5 

3.6 

U10288 

1.53 

0.5^ 

4.4 

0.5^ 


Note.— Col. 1: Name; Col. 2: Flux density C-band (mjy); Col. 3: Uncertainty of flux density C-band, 
in most cases 2% calibration errors; Col. 4: Flux density L-band (mJy); Col. 5: Uncertainty of flux density 
L-band, in most cases 2% calibration errors. 

^The error between the two weightings is larger than the 2% calibration error, and in these cases, we use the 
former. 

"^For UGC 10288, we refer to Paper III. 
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